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Abstract 

We investigate transport properties of a trilayer made of a d-wave supercon- 
ductor connected to two ferromagnetic electrodes. Using Keldysh formalism 
we show that crossed Andreev reflection and elastic cotunneling exist also with 
d-wave superconductors. Their properties are controlled by the existence of 
zero energy states due to the anisotropy of the d-wave pair potential. 
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I. INTRODUCTION 



Transport in superconductor-ferromagnetic hybrid systems has received great attention 
in the last years due to the progress in nanotechnology which made possible the fabrication 
and characterization of various heterostructures [1-5] . 

During the last years phase sensitive tests have shown that the order parameter in cuprate 
superconductors is predominantly of a d-wecve symmetry [6-9]. In d-wecve superconductors 
the zero bias conductance peak (ZBCP) observed in the tunneling spectra results from 
zero energy states (ZES) that are formed due to the sign change of the order parameter in 
orthogonal directions in k space. The ZBCP depends on the orientation of the surface and 
does not exist for s-wave superconductors [10-15]. 

In d-wecve superconductor ferromagnet junctions the ZBCP is suppressed by the increase 
of the exchange field of the ferromagnet [16-19]. This can be understood from the fact that 
increasing spin polarization in the ferromagnetic electrode suppresses Andreev reflection and 
therefore suppresses the ZBCP which is due to the fact that the transmitted quasiparticles 
are subject to the sign change of the order parameter. Moreover in ferromagnet /d- wave su- 
per conductor/ferromagnet double junctions, the quasiparticle current is enhanced compared 
to the normal state because of ZES [20]. Very recently a circuit theory of unconventional 
superconductors has been presented [21]. 

Keldysh formalism has been applied to normal metal-s-wave superconductor junctions 
[22], and in multiterminal configurations where one s-wave superconductor is connected to 
several ferromagnetic electrodes [23-26]. The conductance of multiterminal hybrid structures 
is due to two types of processes: (i) crossed Andreev reflection in which Cooper pairs are 
extracted from the superconductor. The spin-up electron of the Cooper pair tunnels in a 
spin-up ferromagnet and the spin-down electron tunnels in a spin-down ferromagnet; (ii) 
elastic cotunneling in which a spin —a electron from one electrode is transferred as a spin-cr 
electron in another electrode. 

The purpose of the present work is to investigate transport properties of a 
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ferromagnet-d-wave superconductor-ferromagnet double junction via Keldysh formalism. 
We find that crossed Andreev reflection and elastic cotunneling are influenced by the d-wave 
symmetry of the order parameter in the sense that both processes are mediated by zero 
energy states formed for certain orientation of the d-wave order parameter. 

The article is organized as follows. In Sec. II we introduce surface Green's functions. In 
Sees. Ill, IV we describe transport theory and present the results. Concluding remarks are 
given in the last section. 



II. SURFACE GREEN'S FUNCTIONS 

The quasiparticle properties of d-wave superconductors are influenced by interfaces and 
surfaces: due to the anisotropy of the order parameter the quasiparticles that are reflected 
from the surface or transmitted through the interface are subject to the sign change of the 
order parameter. Therefore surface properties are different from bulk properties and we use 
in transport theory the surface Green's functions that take into account the contributions of 
all waves that propagate close to the surface [28,29]. Green's function techniques have been 
used to calculate the conductance of d-wave superconductors near impurities [30,31]. Also 
surface quasiclassical Green's functions have been used in the calculation of the Josephson 
current between d-wave superconductors [32-35] . In our case this Green's function is inserted 
in a three node balistic circuit that is used to describe the transport properties of the hybrid 
structure containing a d-wave superconductor. 

The form of the local retarded surface Green's function matrix for a d-wave supercon- 
ductor is the following, for smooth interface, with momentum conservation in the plane of 
the interface 



9 f 



(1) 



y/ 9 J 

It obeys the Eilenberger equation [36] and satisfies the normalization condition = 1. It 
can be parametrized as follows 



^ 1 + ab'-' 1 + ab'-' 1 + ab' 
where the a{x, 6) and b{x, 6) satisfy the Riccati equations [37] 

hvF cos(e)^ - 2iEa + A*a^ - A = (3) 
ax 

hvF cos(9) ^ + 2iEb - MP' + A* = 0, (4) 
ax 

where I'i^ is the Fermi velocity. We assume for simphcity that the gap function A is constant. 
Then the spatially independent a, b functions are found as 

a{0,e)^ A^f) 

i{E - e.sgnE) 

b{o,e) = , (6) 

where 9 is the angle between the normal to the interface and the trajectory of the quasiparti- 
cle. A_|_(6') = A{9) (A_(6') = A(7r — ^)) is the pair potential experienced by the quasiparticle 



along the trajectory 9{t:—9) and e± = \J ~ ■^±(^)^- In case of dj;2_j^2-wave superconductor 

A(e)=AoCos[2(^-/3)], (7) 

where (5 denotes the angle between the normal to the interface and the x-axis of the crystal. 
Then g and / Green functions are calculated from Eqs. 2. 

In fact in the transport equations we should average over the Fermi surface in order to 
include the details of the order parameter symmetry. We have calculated the density of 
states averaged over the Fermi surface p^^ Reg^^'^{E,9, ) >. In Fig. 2 the density of 
states p^^ is plotted for different orientations of the order parameter /5 = 0, and /5 = 7r/4. A 
ZEP is formed for /3 = 7r/4 due to sign change of the pair potential. A small imaginary or 
effectively dissipative term [5 — 0.01) was added in the energy in order to make this peak 
visible. 
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In order to study the effect of the angular dependence of the transmission coefficient we 
calculate the Fermi surface averaged density of states defined as 

xa;,D^ deD(e)(Reg'''''^(E,e)), 

^ J--K/2 

where D{9) — sin^(^) [38]. We see that the density of states is suppressed when this coef- 
ficient is included in the calculation (see Fig. 2). However we do not expect the transport 
properties to change qualitatively compared to the case where the transmission is indepen- 
dent on the angle (see Fig. 2). In the following we use D{9) — 1. Transport can probe the 
symmetry of the d-wave order parameter if we consider the orientation of the d-wave order 
parameter /3 as a variable. 

The ferromagnetic electrodes are described by the Green's function 

'1.1 

(8) 

p2,2 

where pi^i and p2,2 are respectively the spin-up and spin-down densities of states. 



III. TRANSPORT THEORY 

We use a Green's functions method to describe transport in a system made of two 
ferromagnetic electrodes connected to a d-wave superconductor (see Fig. 1). We first solve 
the Dyson equation which in a 2 x 2 Nambu representation has the following form for the 
advanced (G^) and retarded (G^) Green's functions [39,40] 

qR,A ^ ^R,A ^ ^gR,A ^ |. ^ ^R,A_ (9) 

S is the self energy that contains the coupling of the tunnel Hamiltonian. The tunnel 
Hamiltonian associated to Fig. 1 takes the form 

^ ^ [^^a,a;CQ Cx "I" tx^a^x C-a ~l~ ^6,a;Cj, Cx "I" tx^bCx Cb ■ 
a 

g in Eq. (9) is the Green's functions of the disconnected system (i.e., with S = 0). The 
symbol ® includes a summation over the nodes of the network and a convolution over time 
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arguments. Since we consider stationary transport this conclusion is transformed into a 
product by Fourier transform. G is the Green's functions of the connected system (i.e., with 
E 7^ 0). The Keldysh component is given by [40] 



I + G 



R 



/ + S (g) 



(10) 



The current is related to the Keldysh Green's function [40] by the relation 

e 



(11) 



At this stage of the calculation no explicit angular form of the tunneling matrix elements 
was assumed. The effect of the angular dependence of the transmission coefficient in the 
transport properties was already discussed in the previous section. The elements of the 
differential conductance matrix that we want to calculate are given by 



dl 



(12) 



The principle of the calculation of Qai,aj{Va, H) is similar to the s-wave case [23]. Depend- 
ing on the orientation of the magnetizations in the two ferromagnetic electrodes we can 
distinguish the following cases: 

A. Antiparallel magnetizations 



If the two ferromagnetic electrodes have an antiparallel spin orientation we find for the 
elements of the conductance matrix 



I /l™.2|j. \2 a,a x,x 
ya,a = +47r \ta,x\ Pl,lPg' 



1 



X 



X>AJ)R 

Q ■ 1. |2|j. 1 2 a,a 6,6, A 
^1''^\i^a,x\ I ''6,0; I Pl,l9'2,2 



\2b,b,Rx,x,R 



X- 



x,x,A fx,x,A 



j:x,x,A £ 
}J J 



X>AJ)R 

+ 2i7r|ia^a;| 1^6,0; I Pl\l92,2 



1 1^ \2 b,b,R x,x,R 
J- - 1 16,1 1 5*2,2 9 



X 



J)AJ)R 



2b,b,A'x,x,A 



fX,x,HfX,x,H 1 I J. 
J J [J- - \tb,x\ 92,2 9 
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(13) 



(14) 



(15) 



and 

1 

p, |2|j. |2 a,a b,b rx.x.R fX.x.A / 1 c\ 

ya,b = -47r 1^6,^1 p^pj^ Pl,lP2,2/ (16j 

The expression of the determinant is the following: 

(17) 

and a similar expression holds for Va- Pi]i, P2% Pg'^ ^-^e the density of states of electrodes 
a,b and the superconductor respectively. Contrary to the s-wave case Pg'^ is not zero for 
E < Aq and the term (13) contributes also to the quasiparticle current even below the 
superconducting gap. In the s-wave case there are simple relations between the conductance 
matrix elements (for instance Qa,a = Qa,b) which means that the transport is mediated 
only by Cooper pairs [23]. In the d-wave case such relations are no more valid because of 
the quasiparticle tunneling. Also depending on the trajectory angle 9 the different terms 
contribute to the Andreev current and the quasiparticle current. Moreover in the present 
case the propagators f^^^^MR) ^ g^^^^MR) have a d-wave symmetry. In (16) Gab depends on 
jx,x,Rjx,x,A g^j^j therefore the corresponding matrix element is associated to crossed Andreev 
reflections. 



B. Parallel magnetizations 



If the electrodes have a parallel spin orientation, we find 



ya,a = +47r \ta,x\ Pl,lPg' 



1 



X 



\2b,b,Rx,x,R 



Q ■ 1. |2|j. 1 2 a,a b,b,A 
^1'T^\<^a,x\ V^b^x] Pl,l9l,l 



1 - Kxl'gTf''''' 



X>Aj)R 

+ 2in\ta,x\%J^p'l:t9i!'i'' 



X 



J)AJ)R 



gX,X,RgX,X,R 



(18) 



(19) 



(20) 
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and 

r> /l^2|j. |2|j. |2 ^ (i,a b,b x,x,R„x,x,A /'01 ^ 

t/a,6 = -47r 1*6,^1 ' ' 5' ' ' ■ (^Ij 

The determinant X>r is given by 

= 1 - |iMl'^2:2'''5^'^'^ - Kxl^gtli'^'g'"''''''. (22) 

In (21) Gab depends on g^^^^Rg^^^A and therefore the corresponding matrix element is 
associated to cotunnehng processes. The elements Qb^a, Qb,b of the conductance matrix which 
describe transport through electrode h are derived from the corresponding expressions for 
Ga,a-i Ga,b by the Substitution a bfor the parallel alignment. For the antiparallel alignment 
the following set of substitutions should be made: g2 2^^^^ ^ gi'T^^^\ta, 



IV. RESULTS 
A. Antiparallel magnetizations 

We consider the ferromagnet/d-wave superconductor/ferromagnet double junction shown 
in Fig. 1. For the antiparallel alignment of the magnetizations in the two ferromagnetic 
electrodes the conductance depends on the orientation j3 as well as on the transparencies of 
the interfaces ta^^, h,x- For P — n/A (see Fig. 3(a)) the surface Green's function has a pole 
at £^ = and the conductance (both and ^"^) acquires a ZEP. The conductance Q"^ 
above the gap depends only on the density of states and for large energies it has a finite 
value. Q"''' depends only on crossed Andreev reflection processes and is zero above the gap. 

For P — (see Fig. 3(b)) similarly to the s-wave case no ZES are formed at the interface 
and both and Q"''^ take relatively small values. However the line shape of the conductance 
is V and is determined by Pg^. In the s-wave case the line shape of the conductance is U 
and a peak just below the energy gap exists [23]. In this sense the results for s-wave are 
qualitatively different than for (i-wave with /5 = due to the anisotropy of the (i-wave order 
parameter. 
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To summarize transport for antiparallel magnetizations is due to crossed Andreev reflec- 
tion in which a spin-up electron from one electrode is transferred as a spin-down hole in 
the other electrode, and is influenced by ZES that are formed at the interface due to the 
sign change of the order parameter. The enhancement of the quasiparticle current at = 
for /3 = 7r/4 in the d-wave superconductor ferromagnet double junction has also been found 
recently using the scattering approach [20]. 

B. Parallel magnetizations 

The results concerning the ZES are not modifled qualitatively when the orientation of 
the magnetizations is parallel (see Fig. 1). The conductance above the gap is determined 
mainly by p^^ . For /3 = (see Fig. 4(b)) the results are similar to the case of the antiparallel 
alignment. 

To summarize transport for parallel magnetizations is due to elastic cotunneling in which 
an electron from electrode a is transmitted as an electron in electrode 6, and is influenced by 
the ZES that are formed for certain orientation of the (/-wave order parameter. For /3 = 7r/4 
the interface at large values of the barrier strength becomes transparent due to bound states 
formed because of the sign change of the order parameter in orthogonal directions in /c-space. 
This property does not exist for s-wave superconductors. 

V. RELEVANCE TO EXPERIMENT 

Multiterminal superconductor ferromagnet structures can be used to test the speciflc 
physics associated to the symmetry of a c?-wave superconducting order parameter, because 
transport through the ferromagnetic electrodes has a strong directional dependence. The 
orientation of the electrodes can be used to probe the symmetry of the order parameter. 

The line shape of the conductance spectra is V-like which is a flngerprint of d-wave 
systems. This has already been tested in experiments [11]. 
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Moreover we used a theoretical description that is valid not only in the tunnel regime 
but also for large interface transparencies. We have found a ZEP in the conductances Qa,b: 
Qa,a in the tunnel regime in the two cases of parallel and antiparallel spin orientations, for 
the P — tt/A orientation. 

VI. CONCLUSIONS 

Using a Keldysh formalism we have shown that in the ferromagnet/d-wave superconduc- 
tor/ferromagnet junction, transport is due to crossed Andreev reflection and elastic cotun- 
nehng and is mediated by ZES that are formed at the interface due to the sign change of 
the order parameter. 

We have used the local surface Green's function given by Eq. (1) and we find no particular 
relation between the conductances in the parallel and antiparallel alignments. In the s-wave 
case and for extended contacts it is possible to show that the average current due to crossed 
Andreev reflection in the antiparallel alignment is equal to the average current due to elastic 
cotunneling in the parallel alignment [23,26]. Discussing multichannel effects for d-wecve 
superconductors is left as an important open question. 
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FIG. 1. The geometry of the junction involving the d-wave superconductor and two ferromag- 
netic electrodes of full spin polarization as indicated in the figure by up and down arrows. The 
orientation of the magnetization of the ferromagnetic electrodes can be parallel or antiparallel. 
The excitation at x gives rise to several outgoing trajectories 6. In the figure only one of these 
trajectories is presented. f3 is the orientation of the d-wave order parameter shown also in the 
figure with respect to the direction x. The electron like quasiparticle 2 is reflected as a hole like 
quasiparticle 1 and an electron like quasiparticle 3. The labels a, b in the figure correspond to the 
electrods a, b respectively. 
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FIG. 2. (a) The solid line represents averaged over the Fermi surface density of states for 
P = 7r/4. A well defined ZEP exists. We have put a small imaginary part 5 = 0.01 in the energy uj 
of the Green's function. The dashed line represents the averaged over the Fermi surface with 
an angular depended transmission coefficient, (b) The same as in (a) but for /? = 



14 



Q) 

E 

0) 
0) 
.X 

fS 
E 

Q) 
O 

o 

T3 

o 
O 



0.02 
0.01 


-0.01 

-0.02 
0.03 

0.025 
0.02 

0.015 
0.01 

0.005 


-0.005 



(a) 







(b) 



0.5 



1 
E 



G 



ab 



1.5 



FIG. 3. Conductance Q°'^,Q'^°' for the antiparallel spin orientation of the ferromagnetic elec- 
trodes, as a function of E (in units of Aq) for different orientations of the d-wave order parameter 
(a) /J = 7r/4, (b) /3 = 0. The hopping element is 0.01. 
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FIG. 4. Conductance for the parallel spin orientation of the ferromagnetic electrodes, 

a function of E (in units of Aq) for different orientations of the d-wave order parameter (a) 
= 7r/4, (b) (3 = 0. The hopping element is 0.01. 
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